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DECAY OF FAR FLOW FIELD IN TRAILING VORTICES 


B. S. Baldwin, N. A. Chigier,* and Y. S. Sheaffer 
Ames Research Center 


SUMMARY 


A finite-difference machine code is used in the wake-vortex problem in the quasi-cylindrical 
boundary-layer approximation. A turbulent-energy model containing new features is developed that 
accounts for the major effects disclosed by more advanced models in which the parameters are not 
yet established. Several puzzles that arose in previous theoretical investigations of wake vortices are 
resolved. 


INTRODUCTION 


Methods for reduction of velocities in trailing vortices of large aircraft are of current interest 
for the purpose of shortening the waiting time between landings at central airports. Procedures for 
estimating the response of the vortices to stimuli such as air injection or special flap deflections are 
required. The favorable effect of a special ogee extension at the wing tip has been investigated in 
wind tunnel experiments by Rorke and Moffitt (ref. 1). It is necessary to determine the effects of 
this proposal and possibly several others in the far flow field as well as in the near field that can be 
investigated in wind tunnels. Computer simulation of this problem may eventually produce compet- 
itive flexibility. 

At the present time definitive simulation of the rate of decay and growth in core size due to 
turbulent viscosity has not been achieved. A useful review of previous work on this problem and 
contributions to its solution have been given by Govindaraju and Saffman (ref. 2). Several attempts 
have been made to simulate the effects of turbulence in wake vortices using integral methods, 
notably by Kuhn and Nielsen (ref. 3) and Fernandez and Lubard (ref. 4). Additional useful informa- 
tion resulted from these efforts. The wind-tunnel measurements of Hoffmann and Joubert (ref. 5) 
and Chigier and Corsiglia (ref. 6) have been used extensively and have contributed to the present 
knowledge of turbulent wake vortices. A preliminary calculation of the flow in an isolated vortex 
has been made by Donaldson (ref. 7) based on a sophisticated turbulence theory. The method of 
invariant modeling employed should be useful in the development of future turbulence theories. 

Investigators at the Imperial College in London (e.g., ref. 8) have developed turbulence models 
in which the length scale is governed by an additional differential equation in lieu of assignment of a 
length scale as in Donaldson’s model. Application of this concept to turbulent wake vortices in the 
present paper leads to a significant departure from Donaldson’s prediction of the rate of decay of 
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such vortices. An effect of suppression of turbulence due to curvature of the flow contained in 
Donaldson’s model, but not in that of Launder and others, also plays an important role. A method 
for including these two effects in an extension of PrandtLs turbulent energy model (ref. 9) is 
developed in the present paper. The relatively slow rate of decay of turbulent wake vortices that has 
been observed in flight experiments (ref. 10) is thereby explained. The approach to self-similar 
solutions predicted in earlier investigations by Donaldson (ref. 7) and Baldwin and others (ref. 1 1 ) 
does not occur when both of the above effects are considered. 


THEORY 


Quasi-Cylindrical Flow 


The basic equations in the boundary-layer approximation are taken from Lilley and Chigier 
(ref. 12). Attention has been confined to incompressible flow in this paper and the equations 
reduced to a nondimensional form through division of the velocities by the free-stream velocity 
division of the pressure by p^W^ 2 , and multiplication of the axial and radial coordinates by the 
factor The resulting equations in a nondimensional form suitable for finite difference 

solution are 
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The turbulent viscosity coefficient vj is determined from a turbulence model to be discussed. The 
boundary conditions are 


u(r 2 ,z) = 0 (8) 

v(r 2 ,z) = \I2 w 2 (9) 


vv(r 2 ,z) = 1 


p(r 2 ,z) = 0 

where r 2 is a radius large compared to the size of the core (/-, ). At the center of the core 

u(0,z ) = 0 

v(0,z) - 0 

The initial conditions are 

u(r,z 0 ) = 0 


Hr,z 0 ) = f v (r/r , ) 
w(r,z 0 ) = 1 + f w (r/r t ) 

where r, is the radius at which the circumferential velocity is a maximum. 
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Turbulence Model 

Since detailed experimental information on the flow in turbulent vortices is not yet available, 
methods that depend on evaluation of parameters from the particular experiment at hand cannot be 
used. Instead, a complete model is needed that bears an approximate relationship to the exact flow 
in a manner analogous, for example, to the relationship between simple kinetic theory and the 
rigorous kinetic theory of gases. It has been shown by Launder and others (ref. 8) that Prandtl’s 
1 945 turbulent energy model (ref. 9) works reasonably well in axisymmetric wakes. The basic 
equation can be written (in the notation of A. Townsend, ref. 13) 
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where F is an additional term developed in this paper. The kinematic-turbulent-viscosity coefficient 
vj and local shear are given by 
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For later discussion it is worthwhile to consider the production term of equation (16) in more 
detail. Since 
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the production term can be written 
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In axisymmetric wakes, the mixing length L is assigned a uniform value that is taken to be a 
constant fraction of the width of the wake. However, Launder and others (ref. 8) have developed an 
additional equation from which the mixing length can be computed rather than assigning its value 
(again in Townsend’s notation) 
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Launder and others (ref. 8) and Hanjalic and Launder (ref. 14) have shown that computation of the 
mixing length in this manner leads to improvement of comparisons with measurements in a variety 
of turbulent free shear flows as well as boundary layers. 


The foregoing turbulence model does not apply to turbulent vortices because it does not 
account for suppression of turbulence due to curvature of the mean How that has been disclosed by 
more advanced models such as that of Donaldson (ref. 7). This effect is illustrated in figure 1 , which 
represents a rotating flow in cross section. The Rayleigh stability criterion indicates that if the 
product VR decreases in the outward direction, the flow is unstable against disturbances that would 
interchange fluid between inner and outer regions. Conversely, if the product VR increases in the 
outward direction, there is a damping effect on turbulent eddies produced by the shearing motion. 
This concept can be made quantitative by finding the energy per unit mass of fluid required to 
interchange fluid between the inner and outer annuli, while retaining an unchanged angular momen- 
tum of the fluid transferred 
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The rate of transfer of fluid by an eddy of size L and velocity q 2 is approximated by 
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These relations can be combined to obtain an additional term to be included in the turbulent energy 
equation (18). The new turbulence suppression term is 

„ _ „V_ 8 (VR) , /Tr 

F — 2 g 2 L v *7 ^3 (23) 

This relation applies more generally to flows with curvature if R is taken to be the local radius of 
curvature of streamlines. Present information does not permit determination of a better value of the 
magnitude of this term than that resulting from the above derivation (C 3 = 1 .0). 

It is interesting to note that the suppression term is of the same form and same sign as a term 
arising from the shear that is already present in Prandtl’s energy equation, namely, the last term in 
equation (20). One might wonder whether the curvature suppression effect discussed above is 
already included in Prandtl’s model. However, as can be seen in equation (18), the shear contribu- 
tion to the production of turbulence must always be positive because it consists of a sum of squared 
quantities. It seems clear that Rayleigh instability should enhance the growth of turbulence and an 
opposite stabilizing condition of the mean flow should be capable of stopping the growth entirely, 
as in the suppression of turbulence due to buoyancy forces (e.g., ref. 15). We therefore conclude 
that the term given in equation (23) should be included again even if it does already appear in the 
production term due to a different mechanism (shearing motion rather than curvature of the mean 
flow). 

It may be possible to simulate the suppression of turbulence by making the length L smaller so 
that the effect of turbulent mixing is reduced and the rate of production due to shear decreased. An 
improvement in the agreement of calculations with measurements in curved boundary layers was 
achieved in this manner by Bradshaw (ref. 16). In near-equilibrium situations it may be possible to 
relate that method to the present one. However, the mechanism by which such an effect would 
enter the dissipation equation (eq. (19)) is not clear to the present authors, whereas the origin of the 
energy suppression term (eq. (23)) seems relatively obvious. 

As mentioned earlier, the turbulence model of Donaldson (ref. 7) contains a curvature suppres- 
sion effect that enters through the tensor formalism used in the method of invariant modeling. 
Ultimately, such an approach may be more satisfactory than that employed here. We have been 
unable to relate our method to that of Donaldson and have not been convinced that his treatment 
of the effect is yet in final form. For the present, we prefer simple physical reasoning applied to the 
special type of How under consideration in which the origin of the effect seems understandable. 

The nondimensional forms of the energy and dissipation equations used in the finite-difference 
calculations are 
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The dimensionless local shear .v and length (! are given by 
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The numerical values of the constants (taken from Launder et al . , 
notation. The dimensionless turbulent viscosity coefficient is 

ref. 8 ) are given in the list of 

(^ 7 /^ 00 ) — 

(30) 


Computational Proeedure 

A machine code was developed for finite-difference solution of the foregoing relations accord- 
ing to the explicit predictor-corrector method of MacCormack (ref. 17). Equations (I), (2), (24), 
and (25) are stepped forward in z to obtain values of p w , v, k, and £. At each step, equation (3) is 
integrated by quadrature to optain p. Equation (4) can then be used to evaluate vv. The derivative 
bw/bz is computed by finite-difference and equation (5) is integrated by quadrature to obtain u. In 
the explicit method used to step forward in z, all auxiliary quantities are evaluated in terms of the 
independent and primary dependent variables by means of the remaining formulas. One hundred 
equally spaced mesh points are used for the radial coordinate r. The step size of the axial coordinate 
z is dictated by the requirements of numerical stability. Instabilities arising from the last term in 
equation (1) are dominant according to an analysis based on the methods of Lomax and Bailey 
(ref. 1 8). A workable stability criterion was found to be 

Az < 0.3(Ar ) 2 // max (3 0 

RESULTS 


For isotropic shear flows, the simplified turbulence model described accounts for the major 
effects disclosed by more advanced models in which the parameters are not yet established. It 
therefore applies to turbulent vortices with small axial velocities and can be expected to lead to 
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realistic predictions of changes taking place in the far How field of wake vortices. We have made 
finite difference calculations based on this model in the quasi-cylindrical approximation. Since 
information on the initial values of the axial and radial velocities is unavailable, these velocities have 
been set equal to zero. Thus calculations in this report correspond to an isolated vortex extending 
from — 00 to +°° through the transformation t = Z/W^. 

In the middle of figure 2 are shown computed circulation profiles with free-stream velocity 
Woo, circulation at large radius T^, and wing chord c assigned values corresponding to the Cherokee 
flight data of McCormick, Tangier, and Sherrieb (ref. 10). The calculations were started at z = 0 
with an initial circulation profile corresponding to that of a self-similar laminar vortex (Lamb 
vortex). As shown in the profile at 34.5 chord lengths behind the wing, an overshoot develops in the 
circulation and moves outward. In the inner regions where the circulation rises steeply, the curva- 
ture suppression results in reduced mixing so that changes take place more slowly than in the outer 
regions. The corresponding velocity profiles appear at the top of figure 2. 

At the bottom of figure 2, a typical turbulent intensity profile at a downstream station 34.5 
chord lengths behind the wing is shown. The solid curve represents the computed turbulent inten- 
sity and the dashed curve indicates the equilibrium values that would be present if the production 
and dissipation terms in the turbulent energy equation (eq. (18)) were balanced to zero. The 
equilibrium values are large at a radius of 40 in. (1.016 m) because the gradient of circulation is 
negative in that region as can be seen in the middle curves of figure 2. This negative gradient leads to 
an enhancement rather than suppression of turbulence as indicated by equations (18) and (23). 
However, there is a lag of the actual turbulence below the equilibrium values. The turbulence is 
weak inside the radius of maximum velocity so that relatively slow changes in velocity occur there. 
The local peak in turbulent intensity near the radius of maximum velocity can be attributed to 
changes outside this radius which enhance the local shear leading to production of turbulence. 

Figure 3 contains profiles of the computed mixing length at three axial stations. The calcula- 
tion was started with a uniform mixing length at z = 0. However, the values quickly change within 
10 chord lengths and become highly nonuniform across the profile. In the region of circulation 
overshoot where the turbulence is enhanced, the eddies become quite large. The eddy size quickly 
decreases to very small values inside the radius of maximum velocity where the turbulent intensity 
is also small due to curvature suppression. In figure 4 are shown turbulent viscosity profiles. These 
results indicate that there is considerable turbulent mixing in the outer part of the vortex, but very 
little near the point of maximum velocity. 

Figure 5 contains estimates of the decay of maximum circumferential velocity according to 
several turbulence models. The circles represent the flight data of McCormick and others (ref. 10). 
The lower curve indicates a typical result from previous investigations (refs. 4 and 1 1 ) in which the 
suppression of turbulence due to curvature of the flow is not taken into account. In the middle 
curve, this effect was included but the mixing length was taken to be a constant fraction of the 
radius of maximum velocity. In that case the calculation indicates appreciable encroachment of 
turbulence into the core by diffusion as in Donaldson’s calculation (ref. 7). However, if the variable 
mixing length is accounted for by integration of equation (21), such diffusion is inhibited by the 
very small mixing length that exists near the radius of maximum velocity. The top curve is based on 
the complete model with all large effects taken into account. 
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There arc several possible explanations for the discrepancy between the upper curve and the 
flight data represented by the circles. Calculations based on different initial circulation profiles may 
lead to better agreement. Also, it is known that the level of ambient atmospheric turbulence affects 
the rate of decay. Rapid decays occur in rough air, possibly due to time dependent distortions of 
the circulation profiles introduced by the buffeting. In calm air, slow decays take place that would 
be in better agreement with calculations not including such effects. The important point is that the 
present simplified turbulence model predicts slow decays, in contrast to previous results that were 
in gross disagreement with flight data. This model indicates that the level of turbulence inside the 
radius of maximum velocity is much lower than previously supposed. The initial roll up is more 
nearly inviscid than many of us believed, particularly in the inner regions. Two basic mechanisms 
are responsible for this phenomenon. One is the suppression of turbulence due to curvature of the 
flow (eq. (23)). The other is the variable mixing length (fig. 3) which becomes small near the radius 
of maximum velocity, thereby inhibiting the inward diffusion of turbulence. These effects lead to 
small turbulent mixing in the inner regions of the vortex, so that the approach to self-similar 
solutions noted in previous investigations (refs. 7 and 1 1 ) does not occur. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif. 94035, October 12, 1973 
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Figure 3.- Mixing length = 132 ft/sec (40.2 m/sec), r oo = 112 ft 2 /sec (10.4 m 2 /sec), 

c = 5.26 ft (1.60 m)). 
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Figure 4.- Turbulent viscosity profiles (K' eo =132 ft/sec (40.2 m/sec), r M = 112 ft 2 /sec 

(10.4 m 2 /sec), c = 5.26 ft (1.60 m)). 
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Figure 5.— Decay of wake vortex (R^ = 132 ft/sec (40.2 m/sec), V ^ = 1 12 ft 2 /sec 

(10.4 m 2 /sec), c = 5.26 ft (1.60 m)). 
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